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ABSTRACT 

We calculate light curves produced by a hot spot of a rapidly rotating neutron star, 
assuming that the spot is perturbed by a core r-mode, which is destabilized by emitting 
gravitational waves. To calculate light curves, we take account of relativistic effects 
such as the Doppler boost due to the rapid rotation and light bending assuming 
the Schwarzschild metric around the neutron star. We assume that the core r-modes 
penetrate to the surface fluid ocean to have sufficiently large amplitudes to disturb the 
spot. For a I' — m core r-mode, the oscillation frequency to rj 2mil/[l'(l' + 1)] defined 
in the co-rotating frame of the star will be detected by a distant observer, where V and 
m are respectively the spherical harmonic degree and the azimuthal wave number of 
the mode, and Q is the spin frequency of the star. In a linear theory of oscillation, using 
a parameter A we parametrize the mode amplitudes such that max (|£e|, |^|) / R = A 
at the surface, where £# and £^ are the 9 and cf> components of the displacement vector 
of the mode and R is the radius of the star. For the V = m = 2 r-mode with u> = 217/3, 
we find that the fractional Fourier amplitudes at u = 2fi/3 in light curves depend on 
the angular distance 9 S of the spot centre measured from the rotation axis and become 
comparable to or even larger than A ~ 0.001 for small values of 9 S . 
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1 INTRODUCTION 

. ■— , Accretion powered millisecond X-ray pulsars in low mass X-ray binaries (LMXBs) show small amplitude, almost sinusoidal 
X-ray time variations, the dominant periods of which are thought to correspond to the spin periods of the neutron stars 
J"j , (e.g., Lamb & Boutloukos 2007). Lamb et al (2009) argued that the millisecond X-ray variations are produced by an X-ray 
emitting hot spot located at a magnetic pole of the rotating neutron star, and that so long as the center of the hot spot is only 
slightly off the rotation axis the X-ray variations produced will have small amplitudes and become almost sinusoidal. They 
also suggested that if the hot spot is located close to the rotation axis, a slight drift of the hot spot away from the rotation 
axis leads to appreciable changes in the amplitudes and phases of the X-ray variations. Lamb et al (2009) pointed out that a 
temporal change in mass accretion rates and hence the radius of the magnetosphere, for example, can cause such a drift of 
the hot spot. 

It is now well known that neutron stars can support various kinds of oscillation modes (e.g., McDermott et al 1988). 
There have been, however, no observational evidences that definitely indicate global oscillations of the stars, probably except 
for quasi-periodic oscillations observed in the tail of giant X-ray flares from Soft Gamma-ray Repeaters (SGRs) (Duncan 
1998; Israel et al 2005; Strohmayer & Watts 2005, 2006; Watts & Strohmayer 2006), which are believed to have a global 
magnetic field as strong as B >10 14 G at the surface (e.g., Woods & Thompson 2006). The QPOs observed in SGRs may 
be caused by damping oscillations excited when a sudden restructuring of a global magnetic field in the neutron star takes 
place. For any global oscillations of a neutron star to become observable, mechanisms are needed that excite the oscillations 
to have amplitudes large enough to produce observable variations in the radiation flux. The e mechanism can be an example 
of such excitation mechanisms for low frequency g-modes and r-modes in the surface fluid layers of mass accreting neutron 
stars (e.g., Strohmayer & Lee 1996; Narayan & Cooper 2007). In fact, the r-modes propagating in the surface fluid layers 
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of mass accreting neutron stars have been proposed for the burst oscillations observed during type I X-ray bursts in low 
mass X-ray binaries (Strohmayer et al 1996, 1997; Heyl 2004, 2005; Lee 2004; Lee & Strohmayer 2005). It is also argued that 
retrograde oscillation modes can be destabilized by emitting gravitational waves if they satisfy the frequency condition given 
by < uj/Q. < m (Friedman & Schutz 1978), and that r-modes, which are retrograde modes, are to become unstable to the 
gravitational radiation reaction (Andersson 1998), where Q is the spin frequency of the star, uj is the oscillation frequency 
observed in the corotating frame, and m is the azimuthal wave number of the mode. .R-modes in the fluid core of rapidly 
rotating neutron stars have been a subject of intensive studies (e.g., Lindblom, Owen, & Morsink 1998; Owen et al. 1998; 
Lockitch & Friedman 1999; Yoshida & Lee 2000ab, 2001) and are now regarded as a possible candidate for oscillation modes 
of a neutron star that are excited to produce observable effects. 

If a hot spot on a rapidly rotating neutron star can produce clean light curves without any strong harmonics of the 
spin frequency, it may be possible to use the light curves as a probe into oscillation modes that are excited to periodically 
disturb the spot so that the periodicities of the modes manifest themselves in the light curves. In this paper, we calculate 
light curves by a hot spot taking account of the disturbances by r-modes in the fluid core, which are assumed to be excited by 
emitting gravitational waves. Although the amplitudes of the r-modes are confined into the core, in the presence of a global 
magnetic field the amplitudes may penetrate to the surface fluid ocean (e.g., Lee 2010). In general, the radial component 
of the displacement vector of the r-modes at the surface have amplitudes much smaller that the horizontal and toroidal 
components, and hence the temperature variations at the surface produced by the r-modes would be too small to have 
observable amplitudes. If the horizontal component of the r-modes is large enough to deform the hot spot appreciably, on the 
other hand, the periodic deformation of the spot may manifest itself in the light curves produced by the hot spot. §2 is for 
method of solution, and numerical results are given in §3, and §4 is for conclusion. 



2 METHOD OF CALCULATION 

We consider a rapidly rotating neutron star having a hot spot on the surface and assume that r-modes in the fluid core are 
destabilized by emitting gravitational waves (Friedman & Morsink 1998; Andersson 1998). For a distant observer, light curves 
produced by the hot spot will then be observed to have periodic flux variations with the dominant period equal to the spin 
period of the star. If the core r-modes excited by emitting gravitational waves can penetrate the solid crust to have sufficient 
amplitudes at the surface to give periodic disturbances to the spot, it is expected that the periodicities due to the core r-modes 
will be contained in the light curves from the hot spot. It is the aim of this section to present a method of calculation of light 
curves produced by a surface hot spot which is disturbed by a core r-mode. 

It is convenient to use two Cartesian coordinate systems (x,y,z) and (x',y',z') and to assume the distant observer is 
in the x-z and x'-z' planes, where the z-axis is the spin axis of the star and z'-axis is pointing to the distant observer and 
y — y' ■ We denote by (r, 9, <f>) and (r', 9' ,<j)') spherical polar coordinates associated with the Cartesian coordinates, where the 
z-axis (z'-axis) corresponds to the axis defined by 9 — (9' = 0). We denote by i the inclination angle between the z-axis and 
z'-axis. If we assume Schwarzschild metric to calculate photon trajectories around the spinning neutron star, for a photon 
emitted from a point (R,9',4>') on the surface of the star and reaching the distant observer, the angle 9' is given by (see 
Pechenick, Ftaclas & Cohen 1983) 



,,./2 
R 



dr' rl 1 {. r^l- 1 / 2 

I v-v 



b 2 r'" 



where b is the impact parameter given by 
R 



sin 5 (2) 



V 1 - r a/ R 

where 5 is the angle between the surface normal vector n and the direction vector I, measured by a non- rotating observer at 
the stellar surface, of the photon that reaches the distant observer, and r g = 2GM/c 2 is the Schwarzschild radius, and R and 
M are the radius and mass of the star, and G and c are the gravitational constant and the light velocity, respectively. Note 
that the three vectors n, I, and k! are coplanar, where k! is the unit vector along the z'-axis. The observed differential flux 
cLFe may be given by (e.g., Poutanen & Gierliriski 2003) 

dF E = I E dO = ^~ rg/i V 4 (*) P4j cos SdS, (3) 

D z 11V ' a cos 9' 

where Ie and dO are the intensity of radiation at energy E and the solid angle seen by a distant observer, D is the distance 
to the observer from the centre of the star, I E (S) is the intensity of radiation at energy E into the direction angle 8 measured 
from the surface normal, and dS is the area of a surface element on the surface and the hatted quantities indicate those 
defined in the frame corotating with the star, and we have used cos SdS = cos SdS (e.g., Ghisellini 1999; Lind & Blandford 
1985). Here, r\ is the Doppler factor given by 
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r?= 7(l-/?cosC)' ^ 

where 

Integrating the flux cLFe by photon energy E measured by a distant observer, we obtain 

dF — — ^ — -ri cosd- — J (5 )dS, (6) 

where we have used 



7= — p = — - sin fl and cose = — sirusinra. (5) 

a, ^ a, r?7 sin 6*' 



E = riy/l-r g /RE, r/cosS = cosS, and I (S) = I & (S)dE. (7) 

Jo 

Assuming j(j) = Jo is homogneous black body radiation independent of 5 and integrating over the surface in the corotating 
frame, we have 

,-, /, ,^tR 2 f^ sin 5 a! sin 5 . ~ » - ,„. 

F=(l-r,/*)Jo^, 8 ^^«nW9d* > (8) 

where (f , 6, 0) are spherical polar coordinates in the corotating frame of the star, and we assume 6 — and 4> = cf> + Qt with 
Q being the spin frequency of the star observed by a distant observer, and t is the coordinate time at infinity. Hereafter, we 
let (f , 0, 4>) and (x, y, z) denote coordinates defined in the corotating frame of the star. 

To take account of the effects of periodic disturbances due to core r- modes on the spot, we first calculate small amplitude 
oscillations of rotating and magnetized neutron stars in Newtonian dynamics, disregarding general relativistic effects on the 
oscillations (Lee 2010). We introduce a dipole magnetic field given by B = /i m V(cos#/r 2 ) whose magnetic axis is assumed to 
align with the spin axis of the star, where \x m is the magnetic dipole moment. Since the magnetic pressure in the deep interior 
is much smaller than the gas pressure for neutron stars with a magnetic field whose strength at the surface is comparable 
to or less than ~ 10 12 G, we treat the fluid core as being non-magnetic (e.g., Lee 2007, 2010). Since we assume the spin axis 
is the magnetic axis, the temporal and azimuthal angular dependence of oscillations can be represented by a single factor 
e 1 ( m< p+ w *) ; where in is the azimuthal wavenumber around the rotation axis and u) = a + mQ, is the oscillation frequency in 
the corotating frame of the star with a being the oscillation frequency in an inertial frame. Since the angular dependence of 
the oscillations in a rotating and magnetized star cannot be represented by a single spherical harmonic function, we expand 
the perturbed quantities in terms of spherical harmonic functions YJ m (6, <j>) with different Is for a given m, considering that 
the axis of rotation coincides with that of the magnetic field. The displacement vector £ is then represented by a finite series 
expansion of length j max as 



Jmax 



(9) 



where e r , eg, and e^ are the orthonormal vectors in the f, 8, and <j> directions, respectively, and 

V H = e e — +e - — , (10) 

96» sm8d<t> 

and lj = \m\ + 2(j — 1) and l'j — lj + 1 for even modes, and lj = \m\ + 2j — 1 and l'j = lj — 1 for odd modes, respectively, 
and j — 1, 2, 3, • • • , j max - Substituting such expansions into the linearized basic equations, we obtain a finite set of coupled 
linear ordinary differential equations for the expansion coefficients Si-, Hi. , and T)/ (e.g., Lee 2007, 2010). When the angular 

dependence of the displacement vector at the surface is represented by the functions Ej (6) defined by 

Jmax 

E r (8) e im * = J2 S h ( R ) Y h (*> 4>) (11) 



j'=i 

3a 



E e (8) e im * = e e ■ £] ]H h {R)VnY t ™{8, $) + T t , (R) e r x V H 1^(0, < 



3 = 1 

Jn 



-ic* • ^ \H h (R)X7 H Yi™(8, $) + T t , (R) e r x V H Y^(8, 



(12) 



(13) 



we can rewrite the displacement vector at the surface as 
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Figure 1. a r , ag, and H^ versus cos 8 for the I' = m = 2 r-modes calculated for a neutron star model composed of a fluid core, a 
solid crust, and a fluid ocean, where the dashed, solid, and dotted curves are for S r , Sg, and 3^, respectively, and the amplitudes are 

normalized by max(|S r |, |3g|, |H^,|). From the left to right panels, fl = £1/ i/GM/Ffi = 0.1, 0.2, and 0.3, respectively. 



i/R = [S r (§) e r + Eg (§) eg + iS (§) e>] expi (m</> + ut) , 

the real part of which is given by 

Re (£) /R = ]E r (9) e r + Eg (§j eg] cos (m<j) + cut) — E^, [Oj e<f, sin (m(j> + cut) 



(14) 



(15) 



Note that m<j} + iot — mfi + at. In Figure 1, the functions ! 



■T, — g 



and Ea, of the I' = m = 2 core r-mode calculated for a 0.5Ma 



neutron star model composed of a fluid core, a solid crust, and a fluid ocean are plotted, from the left to right panels, for 
f2 = Q./ \jGM/R 3 — 0.1, 0.2, and 0.3, where we have assumed Bo = 10 10 G with Bo being the strength of a dipole magnetic 
field at the surface (see Lee 2010), and in each panel the amplitudes of the functions are normalized by max(|H r |, |Eg|, |E^,|). 
As shown by the figure, since the amplitudes of the function E r at the surface are much smaller than those of Eg and E,/, for 
the r-modes, we neglect the term E r in the following. In a linear theory of stellar oscillations, the amplitudes of the oscillations 
are indeterminate, and we have to treat the amplitudes as a parameter. In this paper, using a parameter A we normalize the 
oscillation amplitudes such that max (|He(^)|, |S^($)|j = A in < 8 < ir. 

Let us write as r s — Rn s the vector pointing from the stellar centre to the centre of a circular hot spot in the co-rotating 
frame, where 



h s — sin 9 S cos 4> s i + sin 6 S sin <f> s j + cos 6 3 k, 



(16) 



and 9 S denotes the colatitude of the spot centre measured from the spin axis of the star, and i, j, and k are the orthonormal 
vectors in the x, y, and z directions, respectively. If the angular radius of the circular hot spot is equal to a, we have for 
points r — Rn on the spot 

n s ■ n > cos a, (17) 

where n — sin 9 cos 4>i + sin 6 sin <j>j + cos 9k. The outer boundary of the circular spot is given by h s ■ ht = cos a, that is, 

sin 9b sin 6 S cos \d>i, — 4> s ) + cos 9b cos 9 S = cos a, (18) 

where fib = sin 9b cos cj>bi + sin @b sin (j>bj + cos 9bk. If the hot spot is deformed by an r-mode having the displacement vector 
£, the outer boundary of the spot is approximately given by 



(r d ,8 d ,<p d ) = (R,9 b + S9,<Pb + S<P) 

where 

88 = &(R, 0b, 4>b)/R = Eg(8 b ) cos(m0 6 + cut) 

and 

$<fi = £<j>{R,9b,<l>b)/(Rs'm9 b ) = -E^(9 b ) sm(m4> b + ujt)/sin9b 

where \£/R\ <C 1 is assumed. 

If using n we define a vector n± that is perpendicular to n s as 



n± — h — (n • hs) h s 
we have 



-Ye'e + Xe s d 



(19) 
(20) 
(21) 

(22) 
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Figure 2. Plots of (Xj,, Yj,) and (Xd, Yd) for the spot of a = 20° computed using the functions He and 3^ for the I' = m = 2 r-modc 
at f2 = 0.1, where 9 B = 10°, 20° and 30°, from the left to right panels, and we assume uj = 2f2/3 and A = 0.1. Here, the black curve 
is for the non-perturbed spot (X(,,Y(,) and the blue, red, and green curves are for the perturbed spot (Xd,Y d ) at Qt/(2ir) = 0, 1, 2, 
respectively. 



X — sin 9 sin I 



Y = 



, sin t/ cos i 



sin t/ s cos ( 



(23) 



where e| and ej are the orthonormal vectors in the 9 and </(> directions and are perpendicular to e. s r — n 3 . For n — fib that 
satisfies equation (18), we obtain 

Xl + Y b 2 = 1 - cos 2 a, (24) 

where 



Xt — sin 0;, sin ( 



Yi — — cos 9 S sin # b cos 



+ sin 9 S cos # b 



For h — fid, on the other hand, we have 
X d = sin(4 + <50) sin U b + S</> - $„) , Y d = 



cos # s sin(^b + 89) cos ( <£;, + <50 — S ) + sin 9 S cos(9b + i 



(25) 



(26) 



A plot of (Xd, Yd) in a plane may be regarded as a projection of the outer boundary of the deformed spot onto a plane 
perpendicular to the vector n s . Examples of the plots (Xb, Yb) and (Xd, Yd) are given in Figure 2 for fl = 0.1 and in Figure 3 
for fl = 0.3, where for the spot of a — 20° we have assumed 9 S = 10°, 20°, and 30°, from the left to right panels in each figure, 
and A = 0.1 for the amplitudes of the functions He and H^ of the I' — m — 2 r-modes given in Figure 1. Note that we have 
assumed for simplicity that the oscillation frequency of the l' — m = 2 r-mode is exactly equal to uj = 2mQ/[l'(l' + 1)] = 2S1/3 
in the co-rotating frame. The deviation of the r-mode frequency from the analytic formula 2mQ/[l'(l' + 1)] for neutron 
stars depends on various parameters such as the spin frequency, the solid crust thickness, the equation of state, the thermal 
stratification in the core, the general relativistic effects and so on, and since it is not our main concern here to investigate how 
the light curves depend on the frequency deviation, we have simply assumed ui — 2mQ/[l'(l' + 1)] for the modes. The figures 
show periodic deformation of the spot shape, caused by the I' = m — 2 r-mode, as seen in the co-rotating frame of the star. 
As shown by the middle panels of the figures, for a = 20° and 9 S = 20° the outer boundary of the spot touches the pole of 
9 — 0, at which He and H0 vanish. Since the amplitudes of the functions He and H0 show a sharp increase around the poles 
for Cl = 0.3, the outer boundary of the spot is largely deformed when 9 3 ~ a. 



3 NUMERICAL RESULTS 



0.1 and 

,6, 



Examples of light curves produced by a hot spot of a — 20° are plotted as a function of ilt/(2ir) in Figure 4 for il 
in Figure 5 for Q = 0.3, where 9 B = 10° for panel (a) and 30° for panel (b), and we have assumed M = 1.4Mq, R = 10 D cm 
for the neutron star, and uu — 2S7/3 and A — 0.01 for the I' — m — 2 r-modes. Here, we have also assumed that 7o is constant. 
The figures show that the amplitudes of SF = (F — F m )/F m increase when the angular distance 9 S of the spot centre or the 
inclination angle i increases, where F m is the mean flux. An increase in the spin frequency also tends to increase the amplitudes 
of the variations SF through the Doppler factor particularly for large values of 9 3 and i. For a oscillation amplitude A — 0.01, 
we can clearly see the r-modes cause periodic modulations of the amplitudes of SF. 

Using light curves F(t), we calculate the discrete Fourier transform a,j (j — —N/2, ■ ■ ■ ,N/2 — 1) defined by 



^F(t fc )exp(27ri/ j t fc ): 



(27) 
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Figure 3. Same as Figure 2 but for Q = 0.3. 
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Figure 4. Light curves F/F m produced by a hot spot of a = 20° for Q = 0.1 as functions of Qt/(2w) for 9 S = 10° (panel a) and 30° 
(panel b), where F m is the mean flux, and the dashed, dash-dotted, dotted, and solid lines are for the inclination angle i = 10°, 30°, 50°, 
and 70°, respectively. Here, we assume M = IAMq, R = 10 6 cm for the neutron star, and u> = 2H/3 and A = 0.01 for the U = in = 2 
r-mode. 



where N is the total number of sampling points in the time-span AT, tk — kAT/N, fj — j/AT, and \aj\ = \o—j\ for a 
real function F(t), and k and j are integers. For light curve calculations, we use N = N1N2 with iVi = 2 6 and JV2 = 2 5 so 
that we have the Nyquist frequency um v — /(2<5t) with St — P 3 /Ni and the time-span AT = P S N2 where P s = 27r/il is the 
spin period. In Figure 6 plotted are the fractional Fourier amplitudes aj/ao, which is proportional to the fractional rms, as 
functions of o-j/ft with ex, = 27r/j for light curves calculated assuming M — 1.4Mq, R — 10 6 cm, and Cl = 0.3. The dominant 
peak of aj/ao appears at Oj — Q, due to the spin frequency $1 of the star, which we may call the fundamental, and there also 
appear weaker peaks at aj — 20. (first overtone) and 3f2 (second overtone). Because of the periodic modulations caused by 
the I' — m — 2 r-mode having the frequency ui = 2mQ,/[l'(l' + 1)] = 2S7/3 in the co-rotating frame, we also have a noticeable 
peak at Oj — 2H/3. Note that a distant observer will detect the r-mode frequency measured in the co-rotating frame, instead 
of that in an inertial frame, because we are seeing waves restricted to the spot co-moving with the star. Although the peak at 
Oj — 2Q/3 is almost insensitive to the inclination angle i, the peak at Oj — Q becomes higher as i or 8 3 increases (e.g., Lamb 
et al 2009). It is reasonable that the peak at aj = 2J1/3 is approximately proportional to the amplitude parameter A but the 
peak at Oj — Q is insensitive to it. We also note that much weaker peaks, which seem to be proportional to the parameter A, 
are found at aj = fcfl ± j'2il/3 with k and j being integers as a result of nonlinear couplings between the frequencies kQ, and 
J2Q/3. 
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Figure 5. Same as Figure 4 but for Q = 0.3. 



The fractional amplitudes aj/ao are plotted versus 9 S in Figure 7 for H = 0.1 and in Figure 8 for Cl = 0.3. In general, 
the peaks at Oj — Q. and 20. increase their height with increasing 9 8 . On the other hand, the peak at aj = 2fi/3 shows rather 
complicated behavior with increasing 9 S depending on the parameter a and the functions H#((9) and H^((9), but as 9 S — > its 
height simply increases to become even higher than the peak at aj — 2Q, because in this limit the periodic deformation of the 
spot shape is the only cause for any periodicities. For the amplitude parameter A = 10~ 3 , for example, the amplitude aj/ao 
at (Tj = 2fi/3 stays less than 10 -3 when 9 S >10° for both cases of O — 0.1 and 0.3, but it becomes as large as ~ 0.01 for small 
values of 9 S . We also note that although the peaks at aj = O and 2Q. are dependent on the angle i, the peak at oj = 2fi/3 is 
almost insensitive to i because no effects of the velocity field due to the r- modes are included. Figure 9 gives plots of (Xb, Yb) 
and (Xd,Yd) for the case of Q = 0.3, a — 20°, and 9 3 — 50° corresponding to the deep dip found at 9 3 — 50° in the panel 
(a) of Figure 8. We find that differences between the spot deformations in different phases Qt are rather small compared with 
those found in the plots of (Xb, Yb) and (Xd, Yd) in Figure 3, leading to the dip in aj/ao associated with aj = 2J1/3. We have 
also examined the dependence of the amplitudes aj/ao at aj — 20/3 on the compactness parameter r g /R and have found 
that the amplitudes are only weakly dependent on the parameter. 



4 CONCLUSION 

We have calculated light curves produced by a hot spot of a rapidly rotating neutron star, assuming the hot spot is periodically 
disturbed by the horizontal displacement field of the I' — m — 2 core r-mode, which is assumed to be excited by emitting 
gravitational waves. To calculate light curves, we have taken account of relativistic effects such as the Doppler boost due 
to the rapid rotation and light bending assuming the Schwarzschild metric around the star. We have also assumed that the 
oscillation frequency of the I' — m — 2 core r-mode is exactly equal to W = 2fi/3 in the co-rotating frame of the star. It is found 
that a distant observer will detect in the light curves a periodicity due to the I' — m — 2 core r-mode and that the observed 
frequency will be 2fi/3, the frequency of the mode defined in the co-rotating frame of the star, instead of the frequency 4Q/3 
defined in an inertial frame for the mode. The fractional Fourier amplitude aj/ao at aj = 2J1/3 in light curves is approximately 
proportional to the amplitude parameter A which parametrizes the mode amplitude such that max (|Hg($)|, |3^>(0)|J = A for 
< 9 < tv. We find that besides the parameter A the amplitude aj/ao at aj = 2Q/3 depends on the parameter 9 3 , but it is 
almost insensitive to the parameters i and r g /R. The reason for the insensitivity of aj/ao to the parameters may be that no 
effects of the velocity field due to the r-mode on light curves are taken into account of in the calculations and the periodic 
deformation of the spot shape is restricted to a small area on the surface. For A = 0.001 to 0.01, the amplitude aj/ao at 
aj = 2J1/3 will be ~ 0.001 to ~ 0.01 and becomes comparable to or even greater than that of the first overtone at aj — 2fl, 
particularly for small values of 9 S . 

If we write the oscillation frequency of the r-modes observed in the co-rotating frame of the star as 



w/fi: 



K() 



+ K 2 n 2 + o (n 4 ) 



(28) 
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Figure 6. Normalized Fourier amplitudes aj/ag as functions of aj/Q calculated for light curves produced by a hot spot of a = 20° for 
S = 10° (panel a) and 30° (panel b), where Q = 0.3, and the black and red lines are for A = 0.01 and A = 0.001, respectively, and the 
solid and dotted lines are for the inclination angle i = 30° and 10°, respectively. Here, we assume M = 1.4M© and R = 10 6 cm. 
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Figure 7. Normalized Fourier amplitudes aj /arj as functions of the angular distance 8 S of the spot centre from the rotation axis of the 
star for a = 20° (panel a) and a = 40° (panel b), where the black and red lines are for i = 50° and i = 10°, respectively, and the solid, 
dotted, and dashed lines are for the fundamental Uj = Q, the first overtone aj = 2f2, and the I' = m = 2 r-mode aj = 2£2/3, respectively. 
Here, we assume M = 1.4Mq, R = 10 6 cm, and f2 = 0.1 for the neutron star, and A = 0.001 for the mode. 



the coefficient kq for the r-modes with V and m is simply given by 
No = 2m/ [I' (l' + l)] , 



(29) 



and the coefficient ki depends on the physical properties of neutron stars such as the equation of state and the deviation 
from the isentropic stratification in the core (e.g., Yoshida & Lee 2000a,b). Since the neutron star core is nearly isentropic 
such that iV 2 ~ with JV being the Brunt- Vaisala frequency, we only have to consider the I' = m r-modes, for which we have 
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Figure 8. Same as Figure 7 but for f2 = 0.3 




Figure 9. Plots of (Xf,, Yj) and (X<j, Y^) for the spot of o = 20° and 8 S = 50° , where we have used the functions E§ and H^ computed for 
the V = m = 2 r-mode at H = 0.3 and have assumed w = 2S1/3 and A = 0.1 for the plots. Here, the black curve is for the non-perturbed 
spot (Xj,Yj) and the blue, red, and green curves are for the perturbed spot (X4, YJj) for fii/(27r) = 0, 1, 2, respectively. 



oj ~ kqQ, = 2f2/ (771 + 1), and we obtain the frequency u) ~ 20/3 for m — 2 in the co-rotating frame of the star. Although 
no detection of periodicities ui ~ 2J7/3 associated with the V — m = 2 core r-mode has so far been reported, any detection 
of periodicities caused by the I' — m core r-modes in the X-ray millisecond pulsation makes it possible to use the frequency 
deviation given by Aui = ui — kqQ ~ K2Q. 3 to derive information about the equations of state and the thermal stratification 
in the core. 

Employing spcae-time metric numerically computed for rotating neutron stars in stead of the Schwarzschild metric, 
Cadeau et al (2005, 2007) have discussed the frame dragging effects on the light curves produced by a small hot spot on 
the surface of rapidly rotating neutron stars. Cadeau et al (2007) have concluded that in most cases the differences in the 
light curves between the choice of metric are smaller than the differences caused by the spherical or oblate shape of the 
initial emission surface. We think this is also the case for light curves produced by a finite hot spot periodically modified by 
the r-mode. We may also guess that the effects of the rotational deformation of neutron stars on the frequency and surface 
wave-pattern of the r-mode, for example, could be more significant than the frame dragging effects for light curves. We believe 
that at the current stage of investigation the frame dragging may be regarded as one of the effects that become discernible 
only after very accurate light curve determination becomes possible observationally and theoretically. 

To use a surface hot spot of a rapidly rotating neutron star as a probe into core r-modes, the modes must penetrate 
the solid crust to have sufficient amplitudes at the surface. The detectability of periodicities due to the core r-modes in light 
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curves therefore may depend on, for example, the thickness of the solid crust, the property of the surface fluid ocean, and 
the strength of the magnetic field and so on. In this paper, although we have used the functions Eg and H^> computed for 
a low mass, cold neutron star model with a thick solid crust, which makes the frequency spectrum simple, it is desirable 
to use the functions computed for more massive neutron stars as well, which may have a thin solid crust and hence have a 
different oscillation frequency spectrum from that for a cold low mass neutron star. For accretion powered neutron stars, it is 
also desirable to use mass accreting neutron star models with a hot fluid ocean and solid crust, with which we can examine 
how the amplitudes of core r-modes at the surface depends on the physical properties of the fluid ocean and solid crust. The 
existence of a magnetic field is another important factor we have to consider, particularly to determine the wave patterns at 
the surface. Although the effects of a magnetic field on the displacement vector £ of core r-modes have been examined only 
for a dipole field whose axis aligns with the spin axis, we need to examine how an oblique dipole field and a field different 
from a dipole one change the surface wave patterns (e.g., Heng & Spitkovsky 2009). 

Applying a weak nonlinear theory, Arras et al (2003) have estimated the saturation level for the r-mode energy _E r - m odc 
by considering nonlinear transfer of energy to the sea of stellar inertial oscillation modes of rotating stars with negligible 
buoyancy and elastic restoring force and without magnetic field. For the saturation energy in the strong driving limit of 
the r-mode, they obtained an estimate given by E T -m O de/(0.5MR 2 Q. 2 ) ~ 10 -6 (^ sp i n /10 3 Hz) 5 , which may amount to the 
mode amplitude parameter A ~ 10 - (^spin/lO Hz) ' , where u Bp i n denotes the spin frequency of the star. If this amplitude 
estimation is correct, the amplitude parameter A for neutron stars spinning at u Bp i n ~ 500Hz becomes of order of 10 _ , which 
is by one order of magnitude smaller than the value A — 10 -3 used in this paper for presentation and may suggest that the 
detectability of the periodic modulation due to the r-modes is not very high even for small values of 8 3 ■ In this paper we have 
been discussing the possibility of using X-ray light curves from a hot spot as a probe into the core r-modes of neutron stars, 
but it is also conceivable that periodicities due to the r-modes, which may shake the surface magnetic field, are contained in 
the pulses, possibly as drifting sub- or micro-pulses, from millisecond radio pulsars. 
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